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Q ■ Abstract 
O 

. We study the spatial power spectra of Nikolaevskii turbulence in one-dimensional space. First, 

O _ 

2^ \ we show that the energy distribution in wavenumber space is extensive in nature. Then, we 



demonstrate that, when varying a particular parameter, the spectrum becomes qualitatively indis- 
tinguishable from that of Kuramoto-Sivashinsky turbulence. Next, we derive the critical exponents 

Q. 

^ . of turbulent fluctuations. Finally, we argue that in some previous studies, parameter values for 

?H \ which this type of turbulence does not appear were mistakenly considered, and we resolve incon- 

sistencies obtained in previous studies. 
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The spontaneous formation of spatially periodic structure in reaction-diffusion systems 
was predicted by Turing in 1952 l| and experimentally confirmed many years later 
The so-called Turing mechanism is now widely accepted, and we can retrieve many papers by 
searching for the keyword 'Turing pattern', including a large number written this century ^J. 
Recently, we found evidence that the Turing instability in oscillatory systems can also cause 
an initially uniform state to evolve into a state characterized by spatiotemporal chaos instead 
of spatially periodic structure P, 0]. This newly identified type of chemical turbulence is 
exhibited by the equation 

Mix, t) = -di [6 - (1 + dif]^ - (1) 

which was derived from a class of oscillatory reaction-diffusion systems by means of a phase 



reduction technique ^. An equivalent equation was proposed by V. N. Nikolaevskii as a 
model of seismic phenomena [2|. The uniform steady state of Eq. (^, = 0, is unsta- 
ble with respect to finite-wavelength perturbations when the small parameter e is positive. 
However, this instability does not lead to spatially periodic steady states, because the equa- 
tion possesses a Goldstone mode, due to its invariance under transformations of the form 
—>■ + const., and the corresponding marginally stable long- wavelength modes interact 
with the unstable short-wavelength modes. As a consequence, spatially periodic steady 
states do not appear, and instead spatiotemporal chaos is realized supercritically ja 
Spatiotemporal chaos exhibiting a similar onset has been observed experimentally in elec- 
trohydrodynamic convection ("soft-mode turbulence") in homeotropically aligned nematic 
liquidcrystals Q| and numerically in Rayleigh-Benard convection under free-free boundary 
conditions In particular, Eq. has been applied to the study of the former type of 
convective system. It is thus seen that this class of spatiotemporal chaos appears in many 
types of physical systems, and for this reason, studying Eq. ((T)) is important. In this paper, 
we study the statistical properties of the spatiotemporal chaos exhibited by Eq. (P) in one- 
dimensional space with periodic boundary conditions. Also, we argue that in some previous 
works on Nikolaevskii turbulence, values of e that are in fact inappropriate for studying this 
type of turbulence were used. 

Equation (P) has two parameters, the bifurcation parameter e and the system size L. 
First, we derive the L dependence of the spatial power spectrum S{q) = (l^qp), where Vq 
is the spatial Fourier transform of f = 2dxip and () represents a long time average. The 
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FIG. 1: Spatial power density spectrum S{q)/L as a function of the wavenumber q for L = 
2^,2^*^,2^^, and 2^^ with e = 0.02. The fact that these plots fall on a universal curve independent 
of L implies the extensive nature of the energy distribution in wavenumber space. 

quantity S{q)/L is plotted as a function of the wavenumber q for L = 2^, 2^^, 2^^, and 2^^ with 
e = 0.02 in Fig. ^ There it is seen that S{q)/L possesses a universal form independent of L. 
This implies that the energy distribution in wavenumber space is an extensive quantity. In 



Ref. 



12| , the Lyapunov dimension and the Kolmogorov-Sinai entropy are studied for Eq. 



in the cases e = 0.2, 0.5, and it is shown that these too are extensive quantities. However, 
these values of e are too large for Eq. (H)) to exhibit the type of spatiotemporal chaos in 
which we are interested, as we show below. 

Second, we consider the e dependence of the spatial power density spectrum S{q)/L. In 
the following, we consider only the single system size L = 2^, because, as mentioned above, 
S{q)/L is independent of L when L is sufficiently large. The peaks of the spectrum broaden 
and merge when e increases, as shown in Fig. El In particular, when e is larger than about 
0.1, the spectrum is qualitatively indistinguishable from that of the Kuramoto-Sivashinsky 
(KS) equation, dtip{x,t) = —9^(1 + dl)ip — {d^ipY This can be understood as follows. 
The spatiotemporal chaos exhibited by Eq. ^ with a sufficiently small e arises from the 
uniform steady state = 0, owing to the interaction between the weakly stable long- 
wavelength modes and the unstable short-wavelength modes. The band of unstable modes 
has a width in wavenumber space of order e^/^, lying on either side of g = 1. Therefore, 
unless e^/^ ^ 1, the weakly stable and unstable bands of modes cannot be distinguished, 
and thus the situation is effectively the same as that for the KS equation, in which case 
chaos arises through interactions among unstable long- wavelength modes. As a condition 
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FIG. 2: Spatial power density spectrum S{q)/L as a function of the wavenumber 
q for several values of e with L = 2^. From top to bottom, we have e = 
0.4,0.2,0.1,0.08,0.06,0.04,0.02,0.01,0.008,0.006,0.004, and 0.002. The inset displays the spa- 
tial power density spectrum obtained from the well-known Kuramoto-Sivashinsky equation for the 
sake of comparison. 

to insure that the unstable and weakly stable bands of Eq. (P) are sufficiently separated, 
we conjecture that e^^^ must be at least one oreder of magnitude smaller than 1. Hence, 
in order to clearly observe the characteristic Nikolaevskii chaos exhibited this equation, we 
believe that the 0(e) < 0.01 is necessary. This leads us to conclude that the value e = 0.2 



and e = 0.5 used in Ref. 



121 and Ref. 



14| (which employs a wavelet decomposition) are too 



large to observe this type of spatiotemporal chaos and that the power spectrum found in 
those works is actually that of KS spatiotemporal chaos. In fact, it is shown below that 
the exponents of the e scaling for the chaotic fluctuations of Eq. do not converge for 
0(e) > 0.1. 

Noting that the spatiotemporal chaos exhibited by Eq. (0) results from the interaction 
between the long-wavelength modes near g = and the short-wavelength modes in the 
unstable band surrounding g = 1, P. C. Matthews and S. M. Cox derived closed-form 
amplitude equations by hypothesizing that the behavior of the system can be described in 
terms of a quantity v taking the form 



e^/^A(X, T)e*^ + c.c. + e/(X, T), 



(2) 



where A and / represent the slowly varying amplitudes of the two sets of modes, and we 
define X = e^/^x and T = e^t [l^. Also, H. Fujisaka et al. derived amplitude equations 
applicable in higher- dimensional spaces using the form Eq. ^ y^]. The validity of this 



form is supported by numerical results that show a/ (f ^) oc e^/^ for e G [0.01,0.1] [l5^. 
However, based on an analysis of the time series of the spatial Fourier amplitude of turbulent 
fluctuations for Eq. (Q) with the fixed value e = 0.0001, Tribelsky and Tsuboi conjectured the 
scaling a/ ^) oc ^^'^ in Ref.Q. Also, in Ref.Q it is shown that ^(t^ oc for e G [0.1, 1]. 
The discrepancy in the latter case seems to be easily accounted for, as it would appear that 



the result for the exponent reported in Ref. 17[ had not yet converged, because the value of 
e used there is too large. However, the situation is not so clear with regard to the apparent 
inconsistency reported in Ref. j^, because the value of e used there is_certainly sufficiently 



small. Furthermore, we believe that the numerical results of Ref. 15| are insufficient to 
establish the validity of the form given in Eq. ((21) for the following reasons: 1. The results 
were obtained for values of e in a range of only one order, e G [0.01,0.1]. It is quite likely 
that this small range is insufficient to yield a clear result for the power law exponent. 2. 



Studying only the order parameter (f^), we are able to examine the validity of only the 
assumed exponent 3/4 for the amplitude of e*^. Verifying the validity of the other exponents 
requires a different approach. 3. Employing the spatial coordinate X = e^/^x implies the 
assumption that the spatial scale of turbulent fluctuations is very much larger than that 
of the fundamental wave e*^. This scale separation is ensured when e^^^ -C 1. Taking 
this condition as implying that e^/^ can be no greater than 0.1, we obtain the requirement 
0(e) < 0.01 to guarantee sufficient separation of scales. Now, to resolve the inconsistency 
among the results of the previous studies and to examine the validity of Eq. 0, we define 
some new order parameters and examine their e dependence both for smaller values of e and 
over a wider range of values of e than in previous studies. 

First, as shown in Fig. El we find that the results for a/ {v"^) converge for 0.1, where 
we have \/Jv^ oc e^/^. These results indicate that if we wish to study the characteristic 
sptiotemporal chaos exhibited by Eq. (^, we must choose a value of e no greater than 
0.1. Then, as seen in Fig. IH we find Ag oc e^/^, where Aq is defined as the width at half 
maximum for the peak centered at g ~ 1 of the spatial power density spectrum. This is 
evidence that the characteristic spatial scale of turbulent fiuctuations is e""'^. Finally, we 
present in Figs. El and IHl the e dependence of the spatial power densities Si^qo)/ L and S{q^/ L 
for the two characteristic modes go and Qc- These are the wavenumbers nearest to g = and 
g = 1 (explicitly, go = 1 x 27r/L and gc = 81 x 1-n jV) that can be realized in our system of size 
L = 1? with periodic boundary conditions. From the figures it is seen that S'(go) oc e^^^ and 
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FIG. 3: e dependence of \J (w^). The three hnes, included for reference, have slopes of 2/4,3/4, 
and 4/4. 
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FIG. 4: e dependence of Ag. The three lines, included for reference, have slopes of 1/4,2/4, and 
3/4. 

^(gc) oc e^. The former relation cannot be seen as clearly as the latter, because, for larger 
wavelength modes, fewer wavelengths are contained in the finite-size space, and therefore 
there is a larger statistical error in the result, which is obtained by integrating Eq. ^ over 
a finite interval of time. Here, substituting v = e"y4(X, T)e*^ + c.c. + e^f{X,T) into the 
Wiener-Khintchine relation, 

S{q) = L [ {v{x)viO))e-''^^dx, (3) 
Jo 

we obtain 

5(1) = L [ e^''{A{X)A{0))dx, (4) 
Jo 

where we have used {AA) = {AA) = {Af) = {Af) = 0. Now, because the characteristic 
spatial scale of turbulent fluctuations is e~^/^, as shown in Fig. HJ we can reasonably assume 
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FIG. 5: e dependence of S{qo ~ 0). The three lines, included for reference, have slopes of 5/4, 6/4, 
and 7/4. 
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FIG. 6: e dependence of S{qc — 1). The three lines, included for reference, have slopes of 3/4,4/4, 
and 5/4. 

{A(X)A{0)) = exp[— e^/^x]. Then, substituting this into the integrand of Eq. (jH), we obtain 

^(1) = Le2"-i/2, (5) 

Comparing this equation and the result displayed in Fig. IHl we find a = 3/4, which is 
consistent with the result found for the exponent of \/Jv^ obtained from Fig. IHl Similarly, 
we find 

S{0) = Le^^~^/\ (6) 

with (3 = 1. Thus, we arrive at the following conclusions. Our results confirm the validity of 
the form given in Eq. Further, they indicate that the result for the exponent of a/ {v"^) 
given in Ref.^3] is erroneous because the value had not yet converged, as conjectured above. 
Equations El and ini imply that the amplitudes of the Fourier modes with wavenumbers 27r/L 
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(~ 0) and 1 defined in Ref. [2| are 27re^/VVI and j ^fl, respectively. Thus, our results 

indicate that the values of the quantity Ref//v^ in Figs. (3) and (4) in Ref. ^ are of order 1 
and 0.1, respectively, because the parameter values used there are e = 0.0001 and L = 27i/p, 
where p = 3.125 x 10^^. In fact, the figures in Ref.jol support this argument. We believe 
that the reason why P. C. Matthews et al. reported that the form Eq. (jSj) is inconsistent 
with the numerical results in Ref. p] is that they missed the —1/2 appearing in the exponents 
of Eqs. El and IHl which is due to the spatial correlation of turbulent fluctuations. 

In summary, we have found that the spatial power spectrum of Eq. (^J) in wavenumber 
space is an extensive quantity. The spectrum for e > 0(0. 1) is qualitatively indistinguishable 
from that of the Kuramoto-Sivashinsky (KS) equation. We obtained the critical exponents 
of the turbulent fluctuations for Eq. (0), and we found that these exponents converge for 
e^O(O.l). Beyond such a value, because the unstable and weakly stable modes of this 
equation are not well separated, it exhibits spatiotemporal chaos of the KS type, not the 
Nikolaevskii type. Therefore, we conclude that the works presented in Refs. jl2[ 14. [l7|. 
where values of e greater than 0.1 were used, in fact studied KS-type chaos. The numer- 
ical results obtained in this paper are consistent with those given in Ref. and with the 
amplitude equations appearing in Ref.ll^. S. Toh reported that a pulse-distributed model 
reproduces the spatial spectrum of the KS equation We believe that that model is 

applicable also to Eq. ((H), for which pulses become more regulary distributed as e decreases. 
The spectrum of the KS equation possesses a wavy structure for large wavenumbers. We 
believe that for the spectrum of Eq. (|T)), the peaks of this wavy structure, which appear at 
q = 1, 2, 3, ■ ■ ■ , become increasingly sharp as e decreases. 
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